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1 Introduction and Summary 

The understanding of the non-perturbative dynamics of pure Yang-Mills theory is one of 
the outstanding problems of theoretical physics. As was realised a long time ago, one of 
the persistent difficulties hampering any progress is the fact that we are studying a gauge 
theory, and it is clear that the fields that we use to define the theory (in the UV) do not 
create physical states, and are thus not the "correct degrees of freedom" . Instead one should 
switch to some set of gauge invariant variables. Traditionally, this has meant the Wilson loop 
operators, and the consequent introduction of the loop space formalism [ll[2j|3]. Although 
in some heuristic sense, this would seem to bear some relation to the expected appearance 
of a string theory, it is extremely difficult to proceed beyond a few basic stages. We believe 
that although the Wilson loop plays a central role in the theory, as its expectation value is 
a useful order parameter for confinement, it should not be thought of as representing the 
"correct degrees of freedom" . 

Instead, building on a constructiorFlof Karabali, Kim and Nair ^4 , we suggest that the 
theory may be discussed in terms of local gauge invariant variables. The idea of reformu- 
lating Yang-Mills theory in terms of local gauge invariant variables is certainly not new 
and many such proposals exist in the literature E] . The unique advantage of Karabali- 
Kim-Nair formulation however is that it is possible to extend actual computations to the 
point when definite quantitative predictions can be made. For example, in [4 string tension 
in pure YM theory in 2 + 1 dimensions has been computed which argees beautifully with 
numerous lattice simulations [SJ [S] . 

Another central element of |4] is the use of a Schrodinger/Hamiltonian approach. Such 
a formalism [lOl [11] is not often used in quantum field theory because it is fraught with 
regularization issues, and it is somewhat of an art to navigate through them. We believe 
however that Hamiltonian formalism together with gauge unvariant local variables make a 
powerfuU combination which may help to uncover many interesting and non-trivial aspects 
of nonperturbative dynamics of Yang-Mills theories. 

In the present paper we wish to make the next step in the Karabali- Kim-Nair program 
and address the question of determination of glueball mass spectrum in pure 2 + 1 dimen- 
sional Yang-Mills theorjQ In principle, masses of glueball states can be extracted from the 
exponential fall-off of vacuum correlators of various gauge invariant probe operators and 
this is the point of view we take in this paper. This technique however requires sufficient 
knowlege of the vacuum of the theory. We believe that it is useful to approximate the vac- 
uum wave functional as a generalizecPl quasi- Gaussian in a certain variable (appropriate to 
a "correct degree of freedom" describing fluctuations around the vacuum) but with a very 
non-trivial kernel which contains information about an infinite number of physical states. 
It is in this sense that a simple idea like a quasi-Gaussian wave functional, reminiscent of a 
"constituent picture" of glueballs ^13] , is capable of encoding a stringy spectrum. We note 



^At the moment Karabali-Kim-Nair construction applies only to 2 -|- 1 dimensional Yang-Mills theory. 
There is however a close relation between the Karabali-Kim-Nair formalism and that of Bars [6|. This will 
be explored in Ref. [7] in the context of pure Yang-Mills theory in 3 -I- 1 dimensions. 

^A brief summary of our main results can be found in our recent publication |12| . 

^This notion will be made precise below. 



in passing that many of these ideas have analogues in condensed matter physics (e.g., su- 
perfluids, superconductors and quantum Hall fluids |14j). The resulting picture of glueballs 
is more reminiscent of open strings, as opposed to the closed-string picture suggested by 
Wilson loops. 

We want to be absolutely clear that the main outcome of our analysis in this paper is 
a proposal for a new non-trivial expression for the vacuum wave functional. As we show in 
Chapter 3 this wavefunctional can be derived by solving the functional Schrodinger equation 
to quadratic order in local gauge invariant variable 9j|^of Karabali, Kim and Nair and 
under a certain extra assumption on the spectrum of the kinetic energy operator. At the 
moment this latter assumption should be considered as conjectural however and, at least in 
principle, it should be possible to prove (or disprove) it by direct computatioij^ 

It would be difficult to judge the validity of our construction were it not for the existence 
of "experimental" lattice data. Pure Yang-Mills theory on the lattice has been studied 
primarily by Teper and collaborators [8j, and information on a significant number of low 
lying states with a variety of Lorentz quantum numbers is available. Although we have 
attempted in this paper to provide as much detail as possible, there remain a number of 
outstanding conjectural results that require further research, but it should be clear to the 
reader that our analytic results agree extremely well with the lattice data. This agreement 
is non-trivial on many accounts, as we will explain below. 

One topic that we will not touch on in this paper is the role of topology. Clarifying 
this point would certainly be of interest. In particular, a variety of (non-gauge-invariant) 
configurations have, over the years, been proposed as relevant to confinement. We should 
note though that there is a crucial entry for topology in confining theories, that of the 
"compactness" of the configuration space leading to a discrete spectrum (TSl HHJ [T7] . The 
interpretation of the kinetic energy operator as a Laplacian on a compact configuration 
space has been discussed previously in ^E\ H] , and is an essential motivational part of our 
analysis in this paper. 

Finally, in principle our vacuum wave functional is good (as a first approximation to 
the lattice data) for any rank N of the gauge group. Thus the ratios of glueball masses 
to the string tension (string tension being different for different ranks of the gauge group) 
approximate the available lattice data reasonably well for any N . However, as was found 
previously by Karabali-Kim-Nair, and is confirmed by our analysis, the string tension com- 
pares best with lattice simulations at large N, and because of this observation, we expect 
that predictions for the glueball masses, as implied by the vacuum wave functional discussed 
in this paper, work best at large N as well. 

We have organized the paper as follows. In Section [2] we provide a review of the Karabali- 
Nair parameterization of 2 -I- 1 dimensional Yang-Mills theory, emphasizing the aspects that 
are most important for our construction. In particular, we focus on the appearance of 
holomorphic symmetry, which plays a central role in the theory. In a subsection, we also 
provide a discussion of the action of spin, parity and charge conjugation in these variables, 
which is crucial for both the construction of the vacuum wave functional as well as the 
correct identification of physical states. We also give some details of the Hamiltonian in 

*This variable can roughly be thought of as magnetic field B. See equation JlSl below. 
^Preliminary analysis of the issues involved is presented in Appendix A. 



these variables, and how the theory can be regulated in a gauge and holoniorphic invariant 
fashion. In Section l3] we discuss our proposal for the vacuum wave functional, its likely 
validity and some of its physical consequences. We then discuss the Schrodinger equation and 
its solution. As was mentioned above this part of the analysis contains many regularization- 
related subtleties which have not yet been completely clarified. The form of the Schrodinger 
equation that we employ should thus be considered preliminary; however, the central point 
of the present paper is simply that the solution of this Schrodinger equation fits lattice data 
remarkably well. Further clarification of these results, we hope, will be uncovered in future 
work. Four Appendices are also provided as supplementary material for this discussion. In 
Appendix A, we give some detail on the computations required for the Schrodinger equation, 
focussing on the difficult regulator issues. In Appendix B, we compute the (divergent) 
vacuum energy, and in Appendix C, we discuss the general solution of the Schrodinger 
equation. It is in the later discussion that we see a crucial result: the only normalizahle 
solution to the Schrodinger equation corresponds to the confining vacuum. In fact, the 
vacuum wave functional can be thought of as interpolating between perturbative regime in 
the UV and the confining physics of the IR. Appendix D clarifies the appearance of glueballs 
as single particle asymptotic states. In Section 4, we discuss invariant probe operators of 
definite J quantum numbers whose correlation functions can be used to extract the mass 
spectrum of glueball states. Because of the non-trivial quasi-periodic structure of the kernel 
appearing in the vacuum wave-functional, we find an infinite number of such excitations. 
In this section we also investigate more fully the spectrum of low-lying glueball states, their 
approximate degeneracies and the corresponding Regge trajectories. We conclude this paper 
with a series of open questions that we have left for future work. 

2 Formalism 

Consider pure SU{N) Yang-Mills theory in 2 + 1 dimensions in the Hamiltonian gauge 
Aq = 0. It is convenient to parameterize the spatial coordinates in terms of complex 
variables z = Xi — 1x2 and z = xi + 1x2, and we will write d = d/dz, d = d/dz. The spatial 
components of the gauge field may be written 

71=^(^1 + ^2), A^]^{A,^zA2). (1) 

We will frequently use the notation Ai^2 = —it°'Ai2 where t"" are N x N matrices repre- 
senting the SU{N) Lie algebra [i",t^] = if^''^" with normalization Tr(ri'') = i,5"^ When 
appropriate, we will use a specific complex index notation for the components of these 
matrices, t"'-., with i,j ~ 1,...,N. 

The quantization of this theory can be considered within the Hamiltonian formalism. 
Our approach will be based upon a change of variables, whose many details have been 
worked out by Karabali, Kim and Nair |3]. The Karabali-Nair parameterization is 

^^J = - E 5Af,:a(M-i)„j , A,^ = + ^(Aft-i),a9(Af t)^. . (2) 

a a 

Here, M is an invertible complex matrix variable, whose index structure we have denoted 
explicitly; generally, we will simplify notation by dropping explicit indices and ordering the 



expressions appropriately. The tracelessness of the gauge field corresponds to the unimod- 
ularity of M and so M e SL{N,C). Note that a (time independent) gauge transformation 

A ^ gAg-' ~ dgg'\ A ^ gAg-^ - dgg~\ (3) 

where g £ SU{N) becomes simply M i— > gM. Correspondingly, a local gauge invariant 
variable is 

H = ]VPM. (4) 

The definition of M implies a holomorphic invariance 

M{z,z) i-> M{z,z)h^z) 

(5) 
Mt(z,z) ^ h{z)M^z,z) 

where h{z) is an arbitrary unimodular complex matrix whose matrix elements are indepen- 
dent of z. This is distinct from the original gauge transformation, since it acts as right 
multiplication rather than left and is holomorphic. Under the holomorphic transformation, 
the gauge invariant variable H transforms homogeneously 

H{z,z) ^ h{z)H{z,z)h\z). (6) 

The theory written in terms of the gauge invariant H fields will have its own local (holo- 
morphic) invariance. In a sense, we have replaced one local invariance by another. The 
gauge fields, and the Wilson loop variables, know nothing about this extra invariance. We 
will deal with this by requiring that the physical state wave functionals be holomorphically 
invariant. 

One of the most extraordinary properties of this parameterization however is that the 
Jacobian relating the measures on the space of connections C and on the space of gauge 
invariant variables H can be explicitly computed 51 

d^l[C] = a dn[H] e2cASwzw[ff] (7) 

where ca is the quadratic Casimir in the adjoint representation of SU{N) (ca = N), cr is a 
constant determinant factor and 

Swzw[H] = ^ [d'^z Tr (dHdH-^) +-^ fd^x €'"'^Tr{H-^dnHH-^d^HH-^dxH) (8) 
2n J I2n J 

is the Wess-Zumino-Witten (WZW) action for the hermitian matrix field H, which is both 
gauge and holomorphic invariant. The quantum inner product may be written as an overlap 
integral of gauge invariant wave functionals with non-trivial measure 

(1|2) = I dfi[H] e2<=^5"'z"'W ^*^2. (9) 

From many of the above expressions, it is clear that a useful gauge-invariant variable is the 
"current" 

J=^dHH-\ (10) 

TT 



In particular, it is easily established that the standard Wilson loop operator in the funda- 
mental representation of SU{N) may be written in terms of J as 



(11) 



$(C) = Tr Pexp ( (b dz dHR-A = Tr Pexp (— I dz A . 

Under the holomorphic transformations J transforms as a holomorphic connection 

J ^ hJh-^ + — dhh~^ (12) 

while dJ transforms tensorially, dJ i-^ h{z)dJh^^{z). Consequently, dJ plays a central role 
in this theory; indeed, it is closely related to the magnetic field 

B = -—M^-^dJMl (13) 

CA 



As is clear from (12 1, the current is a connection for holomorphic transformations and there 
is a corresponding covariant derivative. For example, for an adjoint tensor field ^qS, we 
have 

[D,4>]^d(t)- — [J,<P] (14) 

CA 

which also transforms homogeneously under (|5|. Furthermore, note that we may define an 
antiholomorphic current 

j=^dHH-^ (15) 

TT 

which is distinct from the adjoint of J; in fact J^ = H^^JH . It is easily shown that 

[D, J] = dJ (16) 

which may be regarded as a reality condition on J . Using this, we see that 

{dJ)^ = dJ^ ^ d{H-\jH) = H-^[D, J]H ^ H-\dJ)H. (17) 

Thus apart from conjugation by H, dJ is essentially real. This is one of the reasons why it 
is possible to write the theory entirely in terms of J. For notational purposes, we introduce 
the covariant Laplacian 

A^«^. (18) 

We note that BD ~ Dd ^ -{t^/ca) dJ. 

We will also need a variety of fields written in the adjoint representation. Suitable 
formulas are such as 

M°^ = 2Tr(rMt^M-i) (19) 

7J"'' == 2Tr{t''Ht''H-^). 

One may check that these expressions do indeed transform under holomorphic transforma- 
tions in the adjoint representation. 

To conclude this section, we summarize by noting that we may parameterize everything 
in this theory, including the path integral measure, in terms of gauge invariant variables, 
but in doing so we encounter a new holomorphic invariance which is not seen by the original 
gauge fields. It is fundamental to this theory that we take this new invariance into account. 
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2.1 Spacetime Quantum Numbers 

We will often be interested in gauge and holoniorphic invariant operators constructed as 
traces of products of dJ and its derivatives. As it will be necessary to classify such operators 
with respect to spacetime quantum numbers J^*^, we pause to discuss these now. 

Here J is the quantum number associated with the spatial 50(2) subgroup of the Lorentz 
group. Clearly, A {A) carries 5*0(2) charge —1 (+1), and derivatives carry spin Jq ~ —1, 
Jq = +1. Upon the KN reparameterization, H will have spin zero, and thus the current J 
carries angular momentum Jj = —1. Consequently, Jqj = 0. For the most part then, the 
spin of an invariant operator will be determined by the net number of derivatives. 

Parity and charge conjugation are determined as follows. By parity, we will mean the 

operation 

p p 

Xi 1-^ Xi, X2 1-^ —X2- (20) 



Thus, under P we have 



P : z ^^ z 

A ^ A 



(21) 



H ^ H- 

A field which is conjugated by H (up to sign) under parity is transforming tensorially 

$ H^ a^H^^^H (22) 



where a$ = ±1. We then find 



and 



[D,[d,n ^ a^H-^[d,[D,nH (23) 



[d,[DM ^ a^H-^[D,[dMH. (24) 



Derivative operators of definite parity will thus be even(odd) linear combinations of these 
two. Taking the sum of the two, we construct 

[AM^\i\DAdM + [dADM]) (25) 

and 

P: [A,$] ^ a$[A,$]. (26) 

Note in particular that A"<I> has the same parity as $. 

Charge conjugation does not act spatially (z ^ z), but we must have 

C : Afj 1-^ ^-^fi ^i] '"^ ^-^fi- (2') 

For M and M ' , a choice of action of charge conjugation consistent with this is 

M,a ^ {M-^)cn Ml. ^ (Mt-i),a (28) 
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which leads to 



H^p ^ {H~')p^ J^-p ^ -J^a. (29) 



C: Kb ^ ^% (30) 



Note that generally if 

then 

C: i[D,4>]U ^ +i[D,cl,^])pa. (31) 

Thus C counts J's (mod 2) that are not inside D's. 

2.2 The KN Hamiltonian 

The standard F A/2+1 Hamiltonian 

Hym ^T+V= / Tr ( g^M-E." + -^B'' ] (32) 



where as usual, the electric fields Ei play the role of momenta conjugate to Ai. The Hamil- 
tonian can be rewritten explicitly in terms of gauge invariant variables. This was worked 
out in detail by Karabali and Nair, and has the collective field form 

nKN[J]^m( f r{x)—^+ f ^''\z,w)—^-4r^\ + ^ I dj^Br (33) 

where 

m = %^ (34) 

ZTT 

and in most cases fi*^^ can be thought of as 

fl^^z, w) = —Di^G{w - z). (35) 

TT 

Here G{w — z) is the ordinary Green's function defined by dwG{w — z) — 5^'^\w — z) and 

D''^^5^''d^~if'^''—J'^{w). (36) 

CA 

In what follows, however, we will consider the action of TLkn on local gauge invariant 
operators, and for such purposes a more general point-split expression for fi"'' will be needed 
(see H] or Appendix A for details). 

The derivation of this Hamiltonian involves carefully regulating certain divergent expres- 
sions in a gauge invariant manner [J] . We note that the scale m is essentially the 't Hooft 
couphng [18]. 

3 Vacuum Wave Functional 

Our basic goal is to determine the masses of some of the lowest lying glueball states. The 
technique that we will use is to consider vacuum correlation functions of invariant operators 
at large spatial separation. Consequently we wish to determine the form of the vacuum 



wave functional. Of course, the determination of the full expression for the ground state is 
an insurmountable problem. Therefore our goal here is more modest: we wish to determine 
only the quasi-Gaussian part of the vacuum wave functional. Two comments about this 
program are in order. 

First, the main outcome of our analysi^in this chapter should be considered as conjec- 
tural. The derivation of the quasi-Gaussian part of the vacuum wave-functional that we are 
going to present below depends in a crucial way on a certain assumptiorrj about the spec- 
trum of the kinetic energy operator. At the moment we cannot prove that this assumption 
is correct, and therefore all results that follow from it should be considered as preliminary. 

Second, the next natural question to ask here is to what extent (or in which regime) this 
quasi-Gaussian expression approximates (if it is an approximation at all!) the true vacuum 
of the theory. Unfortunately, at the moment this question can be answered only a posteriori. 
In particular, as we will see in the next section, estimates of glueball masses based on our 
wavefunction reproduce the lattice data remarkably well. 

In this respect we want to point also that a useful insight on the possible range of 
applicability of our wave functional comes if we take 1 + 1 dimensional Yang-Mills theory 
coupled to adjoint mattei|jas a guide. A version of this theory has been studied in the light- 
cone formulation [19 : in that context, partons play the role of constituents, and numerical 
work on the lowest-lying glueball states showed that the wave functions of the low lying 
" glueball" states have probability very close to one of being an eigenstate of parton number 
operator. The relation between the partons of this theory and the constituent glue in our 
study is not clear, nevertheless if we take the results of [19] as a guide, we then expect that 
the constituent picture which emerges within our scheme (and which is similar to parton 
picture of [12]) should give a reliable description of low lying glueball states. 

After these preliminary discussion we are ready to proceed to details of our proposal. If 
we look at the form of the KN Hamiltonian, we note that if we were to drop the potential 
term (for example, this would happen in the large m limit), then ^ = 1 would satisfy the 
Schrodinger equation. Furthermore, we note that, given the non-trivial form of the measure, 
this wave functional would be normalizable! Consequently, we are motivated to consider 
a ground state wave functional of the form 'i'o — e^, and we would like to determine P. 
In principle, P can be any functional which is gauge and holomorphic invariant, as well as 
invariant under spacetime symmetries. 

Of course, finding an exact expression for P is a very difficult task. However, as we said 
above, we want to find P to quadratic order only. Therefore one might consider a wave 
functional that is Gaussian in dJ 

containing a kernel K that we wish to determine. 



The form of an ansatz (37) is reminiscent of what one usually takes in variational calcu- 



lations. We want to be clear, however, that our approach is not variational in the sense that 



®The mai n o utcome of our analysis is summarized in eqs. (JSSJI and |62l|. 

^See eq. teGk and further comments that follow it. 

^This theory is related to pure glue case in 2 -|- 1 dimensions by dimensional reduction on a circle. 



we are not going to minimize vacuum energy density with this ansatz. Our intension is to 
find the form of the kernel K{dd/rn^) by trying to actually solve the Schrodinger equation 
to quadratic order in dJ. In other words we are attempting to really find an explicit form 
of P to that order in dJ. 



To do this properly, however, requires a further generalization as (37 1 is not invariant 
under holomorphic transformations (pi). To repair this problem, we take the generalized 
quasi-Gaussian 

■K /■;,.../ A 



^o^exp(-^^-—^ldJK[ — ]dJ+...] (38) 



which is explicitly gauge and holomorphic invariant. Here A — {D, d}/2 is the holomor- 



phic covariant Laplacian ( 18 1 mentioned earlier, which depends on J. Note also that the 
argument of the exponential is real and invariant under spin, parity and charge conjugation. 
It contains an as yet arbitrary dependence on A but is quadratic in the commutator dJ. 
From the results of section 2 on discrete symmetries, it makes sense to organize the wave- 



functional in this way. The ellipsis in ( 38 ) would then contain terms of quartic order and 



higher in dJ. As we will see, the Schrodinger equation can also be organized along these 



lines, and it will be sufficient for our purposes to consider the terms explicitly shown in (38 1 
It is convenient to write K{L), with L = A/m^, as a formal infinite power series 



K{L) = Y. c»^" (39) 



n=0 



with as yet unknown coefficients c„. We see that the generalized ansatz (38) corresponds to 
the expansion 

^ = -:r^E^"^ + --- (40) 

n 

of P in terms of gauge and holomorphic invariant operators 

On= fdJ{A")dJ. (41) 

Note that of all possible local operators that might appear in P (and which are formally 



represented by ". . ." in (38) and (40l) we keep only this subset since only these operators 



contain part quadratic in dJ. In other words, even though the ansatz (38) is not Gaussian 
any longer and contains {via A) terms of higher order in J, it certainly does not contain 
all such terms. As was mentioned above, we intend to solve the Schrodinger equation 
to quadratic order in dJ only and thus we keep only those terms which are required for 
consistency. 

It should be noted that we are taking here a basis of local operators]^ As we have seen, 
this basis is natural from the point of view of holomorphic and discrete symmetries. It is 
far from clear whether or not this is the most convenient choice, and in particular, we would 
not expect that this choice diagonalizes the Hamiltonian. We will discuss many of these and 
related issues, particularly in Appendix A. 



'Nevertheless, the form of K that we wiU arrive at later is an infinite power series and is thus non-local. 
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One way to intuitively motivate the expansion in ( 38 1 is to think of dJ as the relevant 



local probes of real physical states. Then there should exist an expansion parameter, which 



is related to the size of the glueballs. The quadratic term in (38 1 should be then interpreted 
as the leading term in the expansion in the inverse of that effective glueball size. This would 
be very reminiscent of the a' expansion in string theory. 

After such preliminary discussion of the wave-functional and before proceeding to the 
Schrodinger equation and its solution, let us say a few words about the (expected) asymptotic 
behavior of the vacuum state. In the UV, the vacuum wave functional should correspond 
to free gluons 



Because of the relation (13 1 between B and dJ, we see that this is Gaussian in dJ. 

From the form of the vacuum wave functional in the UV, it is clear that one cannot 
expect iiT to be a local functional of L = A/m^. However, we will formally write K{L) 



as an infinite power series ( 39 ) and seek a summation which corresponds to a normalizable 
solution of the Schrodinger equation. As we will see the resulting kernel K will have a very 
non-trivial form, which contains much physical information. 

Before proceeding further, let us note that in the small L (IR) limit, we expect the kernel 
to asymptote to a constant value. In fact, we will find that 

,IR , . „„„ ( ^ /" rr„ d2 



As explained in [4 , this wave functional can be thought of as providing a probability measure 
^5^0 equivalent to the partition function of the Euclidean two-dimensional Yang-Mills 
theory with an effective Yang-Mills coupling g^jj = rngyM- Using the results from pO] . 
Karabali, Kim and Nair deduced the area law for the expectation value of the Wilson loop 



operator (11) 



($) ~ exp(-crA) (44) 

with the string tension following from the results of |20] 

N^ - 1 



o" = 9ym " g^ ^ ■ (45) 

This formula agrees nicely with extensive lattice simulations [S], and is consistent with the 
appearance of a mass gap as well as the large N 't Hooft scaling^ 

3.1 The Schrodinger equation 

Now let us return to our discussion of the derivation of the vacuum wave functional. To 
properly discuss the Schrodinger equation, it is necessary to regulate the Hamiltonian, as 

^"It should be noted however that this result is not completely satisfactory at finite A^: the representation 
dependence of this result would not be correct [16II21J : it is however consistent with "Casimir scaling". We 
will not discuss this important issue (as well as matters relating to the center of the gauge group) in this 
paper, but it is one indication that the formulation may only be consistent at large A''. 
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we have discussed briefly above. It is fairly straightforward to see however, what the general 
form of the Schrodinger equation will be. Given our ansatz, we will find 



'Hkn'^o — £^0^0 — 



Eo+ dJ (7^) dJ 



*o- 



(46) 



The (divergent) vacuum energy £"0 can be isolated, and in Appendix B we give its derivation; 
as expected, the leading divergence in the UV is cubic. Next, what we need to do is compute 
the expression that we have labelled by TZ and set it to zero. This will constitute an equation 



for the kernel K. The ellipsis in (46) stands for the terms which are at least cubic in dJ. 



Given our discussion above, we can neglect these terms. 

Let us consider the various terms in TZ. First, the potential term J dJdJ clearly con- 
tributes a fixed constant to TZ. The remainder of TZ will come from the action of the kinetic 



energy operator Tkn ■ 
tary to derive 



and thus 



Given the form of the vacuum wave functional ^0 
,5^0 SP 



S^^o 



5^P 



SJ^{z) 



*r 



it is elemen- 



(47) 



SJ''{z)J''{w) 



5P 



SP 



SJ^{z)J>'{w) 5J^{z) 6J^{w) 



from which we find 



Ti-KN^O 



TKNP + m n''\z,w) 



SP SP 



CA 



SJ'^{z) SJ''{w) TT-m 



*o 



drdr 



*n 



(48) 



(49) 



The second term in brackets is easy to compute. Since we want to solve Schrodinger equation 



to quadratic order in dJ only, it is enough to use ansatz ( 37 1 from which we obtain 



P = 



and 



SP 

JjHz) 



2cATn^ 



CA'm'^ 



dJ K 



as 



dJ- 



dK 



dd_ 



d.r{z) 



Also to this order we have from ( 35 1 



n^'^z^w) 



CA 



S'^^d^Giw -z) + 



(50) 



(51) 



(52) 



and by putting everything together we see that the second term in brackets in ( 49 ) is equal 
to 



CAm 



dJ 



dd 



K' 



rdd\ 



\ vn? I 



dJ 



(53) 



Therefore, we see that contribution of this term into TZ is of the form LK^{L). 

Up to this point our discussion was fairly straightforward and general. The difficult part 
however is the last TknP term in (49). There are several hard issues here, among which 



include the utility of the holomorphic invariant regulator, normal ordering and renormaliza- 
tion issues, and the choice of basis for operators that we have taken (and consequent operator 
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mixing). At the moment we do not have full analytic control of these issuses. Therefore in 
what follows we present some heuristic arguments to motivate our main conjecture in eq. 



(56 1. Further technical details can be found in Appendix A. 

The simplest way to think of the kinetic energy operator T^n (in particular the J5/dJ 
term) is that it acts homogenously on any local operator-valued function of J. According 
to this the action of Tkn on P from Eq. ( 50 ) would simply give us 2mP since P is quadratic 



in dJ. This picture, however, is certainly not complete as can be easily seen by considering, 
for example, a gauge and holomorphic invariant operator 



Oi 



driDOydJ' 



driddyor 



Br 






adb -fd 



rd dJ 



(54) 



Tkn acting on this should give a holomorphic ally invariant. The JS/SJ term will act to 
count the number of J's in each term - this part of Tkn is not holomorphic invariant, and 
thus the rest of Tkn must act so as to restore holomorphic invariance. In particular, for Oi 
an extra quadratic in J term is generated and we expect to obtain 

TknOi^StuOi. (55) 

Motivated by this heuristic argumentrj we expect that the general result should be 

TKNOn^{2 + n)mOn + ... (56) 

where the ellipsis stands for terms of higher order in dJ but the same mass dimension as 
On and which will mix with ©„ under the action of Tkn- As we said above and would like 
to repeat here once again, we do not have complete analytic control of the action Tkn and 



therefore at the moment eq. (56) cannot be directly derived from the KKN Hamiltonian. 



We will discuss these issues at some length in Appendix A. Nevertheless, as will be seen 



from what follows, (56 1 should be considered preliminary and conjectural, but its validity 



can in some sense be justified by the fact that it leads to sensible physical results, which, 
without further apology, we explore below. 



Returning to the derivation of the Schrodinger equation, asssuming (56 1 we can write 

(57) 



1 



now 

TknP= — I cl-J \ o ^ 

It is convenient to write the factor in braces formally as 

Assembling all of these results, we then find 

1 d 



c„{2 + n)L" )dJ. 



7^ = 



CAfn 



2LdL 



[L^K{L)\ + LK"^ 



and if we set TZ to zero, we will finally obtain an equation for the kernel K 



K 



f;;^["-(^>i 



LK^ 



1 =0. 



(58) 



(59) 



(60) 



^^Other possible motivation for eq. \bQ\ comes from comparison with lattice gauge theory. We plan to 
elucidate this issue in a future publication. 
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3.2 The Kernel 



Thus, we have succeeded in reducing the Schrodinger equation to a differential equation ( 60 1 
for the kernel, which is of the Riccati type. We will present complete details of its solution 
in Appendix C, and note the important features here. Although the Riccati equation is 
non-linear, it is easily transformed into a linear second-order equation of the Bessel type, 
and one finds a general solution of the form 

icuMR±W^ (61) 

where C is a constant and J„ (F„) denote the Bessel functions of the first (second) kind. 
As explained in Appendix C, it is remarkable that the only normalizable wave Junctional is 
obtained for C — > oo, which is also the only case that has both the correct UV behavior appro- 
priate to asymptotic freedom, as well as the correct IR behavior appropriate to confinement 
and a mass gap! This solution is of the form 

Vi Ji(4VL) 

This remarkable formula is reminiscent of similar results in related contexts |22| : here, it 
encodes information on the spectrum of the theory, as we show below. We note that this 
kernel has the following asymptotics (where L ~ —fP'/Am'^) 

p^O, K^\; p^oo, K^2ralp (63) 

consistent with confinement and asymptotic freedom, respectivelyFj 



Now using standard Bessel function identities (see Appendix C, Eqs. (184)) we may 
expand 

Ji{u) 4 „ ^ 1 , ,, 

where the 72, n are the ordered zeros of J2{u). For example, the first few zeros I23J of J2{u) 
are J2,i = 5.14, J2,2 = 8.42, J2,3 — 11-62, J2.4 = 14.80, etc. The inverse kernel is thus 

Now, if we regard L ~ ddjrr? , in terms of momentum p we find 

00 -II 

Here 

M„ = ^— . (67) 

'^^Note that the argument of Bessel functions is imaginary; so instead of J„ we have /„ Bessel functions. 
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As we will see in the next section, M„'s can be interpreted as constituents out of which 
glueball masses are constructed. The first few of these M„'s are 

Ml = 2.568m 

M2 = 4.209m 

M3 ^ 5.810m (68) 

M4 = 7.398m 

M5 = 8.980m. 



It is not difficult now to find a Fourier transform of inverse kernel K ^(p). By rewriting 

IS 



2 ^ V P' + M^ 
we immediately obtain 

K-\x ^y) = S^'\x - y) + 2 E ( '^^'^(^ -y^-^ Ko{M,,\x - y\)] (70) 

n—l ^ ^ 

where Kq{x) is the modified Bessel function of the third kind. At asymptotically large 
spacial separations \x — y\ ^ oo this takes the form 

^ 00 

K-\\x - y\) « - ^(M„)ie-^^"l--^l. (71) 

This is a primary result, which will be important in our discussion of correlation functions 
in the next section. It is important to appreciate that one should not think of the kernel 
iiT as a propagator of some propagating field. Rather, as we discuss in Appendix D, it is a 
building block for such propagators. 

Finally, to compute correlators at spatial separation, such as {dJ°'{x) dJ^{y)), we have to 
evaluate a path integral over the H field with the insertion of the non-trivial WZW measure 
as well as | ^§ | . This would be a very difficult task unless we notice that dJ in fact can be 
treated as a free field. Below we present formal arguments supporting this observation; one 
may also consider this question formally based on a study of the path integral measure on 
the space of J's, although we do not present this here. 

We start with the observation that J variable can be obtained from A by complex gauge 
transformation with matrix M^ 

A = -dMM'^ ^ A""^ = M^AM^ -^ - OMUi^ -^ = - — J. (72) 

CA 

Note also that under such transformation 

A = -BM^ -HI^ ^ A*^* = M^^Aff -^ - BM^'M^-^ = (73) 

and we may think of transformation from A, A variables to J variables as "gauge fixing" 
with gauge condition A = 0. This is very similar to standard treatment of 2D Yang-Mills 
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theory [2T in which case after the choice of some axial gauge, say Ai = 0, the action simply 
becomes 

S2D ^ ^ f d^A-, d^A-^ (74) 

and we can treat the only remaining variable A2 as a free field. Given this similarity, it 
should be clear that we can treat J as a free field in the sense that 

{dr{x) dj\y)) = 5'^'^-^K-\x - y). (75) 

Now it is not difficult to derive an expression for the vacuum expectation of the large Wilson 



loop (111 of area A. To leading order 



($(C)) = (TrPexp { — <f j\) -^ iVexp (~^ fd^xd^y {dJ{x) dJ{y)) ) (76) 



or 

o4 f An 



^($(C)) = exp {~'-^^Mi^-^jd^xd^yK-\x y)^ (77) 



and from this expression we see that the leading (5- function in ( 70 1 gives an area law with 



the same string tension (45) as was found before by KKN. As for the rest of terms that 



appear in ( 70 1 , we may notice that 

J d^xU^\x-y)-^Ko{M^\x-y\)] ^ as A ^ 00 (78) 

and therefore these terms will give corrections to the area law behavior which vanish for 
asymptotically large loops. The appearance of the Kq Bessel functions suggests that our 
wave functional secretely knows about some effective abelian vortex configurations which 
model the repulsive two-particle part of the ground state. This is familiar in strongly 
interacting systems in condensed matter physics, where often the ground state of strongly 
coupled systems is modelled as a product of single particle wave functionals (the obvious 
WKB part of our wave functional) times an effectively repulsive two-particle part which 
minimizes the energy of the ground state by keeping the effective quasiparticles (the J- 
"particles") apart. 

4 The glueball spectrum and the Regge trajectories 

The primary purpose of this section is to provide a detailed analysis of analytic computations 
of the glueball mass spectrum and to compare these results with available lattice data [H] . 
We will probe the spectrum of the theory by considering two point functions of appropriate 
invariant operators at large spatial separation. Such operators may be classified by their 
J^ quantum numbers, which we have described in detail in Section ^ For example, 
we will consider the 0++ states which may be probed by the operator Tr (dJdJ) and the 
corresponding correlator is 

{TvidJdJ)^TTidJdJ)y). (79) 
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But before continuing with this analysis, let us make a few comments on this point. 
One expects that the glueballs are extended spatially in some way (indeed, a common 
phenomenological picture is to consider the resonances as vibrational modes of a closed flux 
string). We will derive the masses of these states by probing them with local operators. It 
is likely that this procedure has limited applicability; in particular, masses of higher spin 
states are likely difficult to extract reliably. These limitations also apply to the lattice work. 
Therefore we will confine ourselves to a detailed analysis of spin-0 and spin-2 sectors of the 
theory, and then indicate what the general picture might look like. 

4.1 Spin zero states 

We begin with the O"'"'^ states. As we remarked above, we may probe these states using the 



local operator Tr (dJdJ) . Given our knowledge of the vacuum wavefunctional ( 38 1 we may 



compute the correlator ( 79 1 as 



(Tr(aj5J), TridJdJ)^) ^ {K-\\x ~ yl))' (80) 

This is of course not exact as we have dropped all interactions!^ but we have kept the 
non-local spacetime dependence of the kernel. Together with the expression for the inverse 



kernel (71) this gives 



(Tr (ajaj), Tr(5J9J) J «-— : V (M„M„0'/V(^"+^^'")l^-^l. (81) 

327r|x - y\ ^^^^ 

From the characteristic exponential decay we can now read off the 0++ glueball masses: 

Mo++ = Ml + Ml = 5.14m 

Mo++. = Ml + Afa = 6.78m 

Mo++.. = Ml + M3 = 8.38m (82) 

Mo++... = Ml + Mi = 9.97m 

Afo++.„. = Afi + A'f5 = 11.55m. 

Since m is not a physical scale, we should re- write these results in terms of the string tension. 
Given equation (45 1 presented above, we can immediately write for the lightest 0"*""^ state 



^-5.14,/^^^ (83) 



and similar expressions for the rest of glueballs listed in (82 1. Note, however, that equation 



(45 1 for the string tension carries an explicit dependence on the rank N of the guage group. 
Because of this, M/ ^/a ratios depend explicitly on TV as well. As an example, in Table [l] 
comparison of our predictions for the mass of the lightest 0"'""'" state at different values of 
N {N — 2,3,4,5,6,00) with lattice data is given. From this we see that simple "Casimir 



^^It would be of interest to consider the effect of interactions. For example, one would expect that these 
give rise to decay widths and mixings of states. Such an analysis would also likely clarify the identity and 
role of various expansions that undcrly our construction. 
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Table 1: Dependence of the mass of the hghtest 0++ ghieball on the rank of the gauge group. 
All masses are in units of the square root of the string tension. The percent difference 
between our predictions and lattice data is given in the last column. 



Gauge group 


Lattice 


Our prediction 


Diff, % 


SU{2) 


4.716 ±0.021 


4.732 


0.3 


SU{3) 


4.330 ±0.024 


4.347 


0.4 


SU{A) 


4.235 ±0.025 


4.233 


0.05 


SU{5) 


4.180 ±0.039 


4.182 


0.05 


SU{6) 


4.196 ±0.027 


4.156 


1.0 


SU{oo) 


4.065 ±0.055 


4.098 


0.8 



Table 2: 0++ glueball masses in QCD^. All masses are in units of the square root of the 
string tension. Results of AdS/CFT computations in the supergravity limit are also given 
for comparison. The percent difference between our prediction and lattice data is given in 
the last column. 



State 


Lattice, N ^ oo 


Sugra 


Our prediction 


Diff, % 


0++ 


4.065 ± 0.055 


4.07(input) 


4.098 


0.8 


0++* 


6.18 ±0.13 


7.02 


5.407 


12.5 


0++** 


7.99 ±0.22 


9.92 


6.716 


16 


0++*** 


9.44 ±0.38 14 


12.80 


7.994 


15 


0++**** 


— 


15.67 


9.214 






scaling" encoded in ( 83 1 describes A^-dependence extremelly well even for A^ as low as A^ = 2. 



For this reason in the rest of our discussion of glueball mass spectrum, we will limit ourselves 
to comparison with lattice at A^ = oo only. Similar comparisons at finite A^ should present 
no difficulties. In particular, our results for 0++ states at A^ = oo are given in Table [2] 

Several comments are now in order. First, there are not any adjustable parameters in our 
theory; the ratio of Mq++ to ^y<T is a pure number and is in excellent agreement with lattice 
data! This should be contrasted with the supergravity results also listed in the table for 
comparison and which result from calculations [25] based on the AdS/CFT correspondence 
[27j : in that case, the overall normalization is not predicted but is instead determined by 
fftting to the lattice data, for example, to the mass of the lowest 0++ glueball. Second, note 
that we have been able to predict masses of the 0++ resonances, as well as the lowest lying 
member. Our results for the excited state masses differ at the 10-15% level from the lattice 
results. We note that precisely these masses are more difficult to compute on the lattice 
[28j . and thus the apparent 10 — 15% discrepancy may be illusory. 

Another possible explanation of such discrepancy comes from the observation that lat- 



^*Mass of 0++*** state was computed on the lattice for SU{2) only |25| . The number quoted here was 
obtained by a simple rescaling of SU{2) result. 
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Table 3: Same as Table [2] but with lattice data relabelled as explained in text. 



State 


Lattice, iV ^- oo 


Our prediction 


Diff, % 


0++ 


4.065 ±0.055 


4.098 


0.8 


0++* 


6.18 ±0.13 


5.407 


— 


Q + + ** 


6.18 ±0.13 


6.716 


— 


0++*** 


7.99 ±0.22 


7.994 


0.05 


Q + + **** 


9.44 ±0.38 


9.214 


2.4 



tice values for the masses of 0++** and 0++*** states are in much better agreement with 
our prediction for o++*** and o++**** states respectively (c/. 7.99 ± 0.22 vs. 7.994 and 
9.44 ± 0.38 vs. 9.214.). Also, numerically, lattice prediction for 0++* mass seems to lie in 
between our numbers for 0++* and O"*""*"** states thus indicating that it might be some kind 
of an "average" . These two observations suggest the following picture: the first two excited 
states have not been properly resolved in lattice computations and have been counted as a 
single mass eigenstate "0++*" with an "average" mass of 6.18 ± 0.13\/ct; as a consequence, 
the rest of the excitations should be labelled differently and lattice "0++**" state should be 
labelled as 0++***, while lattice "0++***" state should be labelled as 0++****. Assuming 
this picture to be correct, we give an updated comparison of our 0++ mass predictions with 
lattice data in Table |3] Obviously, the numerical agreement is quite impressive. 

Let us move on to a discussion of the glueball resonances. We may probe these 
states with the operator Tr (dJdJdJ) . We are thus interested in the correlation function 



{Ti-idJdJdJ)^ Tr{dJdJdJ)y) ^ {K~Wx~y\)y 



(84) 



or 



(Tr (dJdJdJ).^ Tr (dJdJdJ) 



y' 64(27r|x-y|)3/2 



J2 {MnM^Mk) 



Z/2 -{KU+M„, + Mk)\x-y\ 



n, 771, fc — 1 



Using this result, we obtain glueball masses which are the sum of three M„'s 



Afo-- =Mi+Mi 
Mq--, =Mi+Mi 
Mo--.. =Mi+Afi 



Ml = 7.70TO 
M2 = 9.34TO 
M3 = 10.95m. 



(85) 



(86) 



These results are compared to lattice and supergravity data in Table |4] We see that our 
predictions are within a few percent of the lattice data, and are much better than the 
supergravity predictions. 

One last comment about correlator in (85 1. Although it has an exponentially falling 



nature, the prefactor is not of the correct form so as to be directly interpreted as a single 
particle pole (for details, see Appendix D.) We believe that this discrepancy can be accounted 
for by normalizing the probe operator appropriately. Note that this issue does not arise for 
the O"*""*" states. 
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Table 4: glucball masses in QCD3. All masses are in units of the square root of the 
string tension. Results of ADS /C FT computations in the supergravity limit are also given 
for comparison. The percent difference between our prediction and lattice data is given in 
the last column. 

State Lattice, N ^ 00 Sugra Our prediction DifT,% 

0~ 5.91 ±0.25 00 05 4 

0~* 7.63 ±0.37 9.34 7.46 2.3 

0—** 8.96 ±0.65 12.37 8.73 2.5 



4.2 Spin 2 states 

As we have seen in the previous subsection, there is a fairly good amount of "experimental" 
lattice data available for spin-0 sector of 2 + 1 dimensional pure Yang-Mills theory. Also 
the quality of those data appears to be good enough thus making direct comparison with 
our predictions possible. Unfortunately this is not the case with higher spin states: in most 
cases only the mass of the lowest resonance with given J^'~^ quantum numbers is known (at 
best) and this, combined with much larger "experimental" (lattice) error bars, makes for 
a much more problematic comparison. The only exceptions are spin-2 states and it is our 
intention to discuss them next. However, given the comments presented above regarding 
the amount and quality of the lattice data available at the moment of the writing of this 
paper, one should not really expect such an excellent agreement as we had in spin-0 sector 
(although, as will be seen below, the agreement is still quite impressive) . 

But before we proceed any further let us remind the reader about one peculiar property 
of 2 + 1 dimensional Yang- Mills theory known as parity doubling. Consider a state |P, J; E) 
of certain parity P, angular momentum J and energy E. Note however that since parity 
and angular momentum operators do not commute (they rather anticommute, {J,P} = 0, 
meaning that parity flips sign under the action of J) this state is not an eigenstate of J but 
of J^ only. Therefore by acting on this state with J we obtain an essentially different state 
(or zero if J = 0) 

\-P,J;E) = J\P,J;E) (87) 

of the same energy E, angular momentum J but opposite parity. This means in particular 
that masses of glueballs with the same J and C quantum numbers but opposite parities 
should be the same. Of course, this argument does not apply to J = states and there is 
no reason to expect parity doubling in the spin zero sector of the theory. 

This phenomenon of parity doubling simplifies our task since all we have to do is to find 
masses of 2++ and 2 resonances; then these results should apply, respectively, to 2 ^ 
and 2^ states as well. For 2++ states we may take the simplest possible spin-2 operator 
namely Tr (d^Jd'^J) and compute the correlator 

(Tr(aVaV)^ Tr(aVaV)^) ^ {d,dyK-\\x-y\))' 
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Table 5: 2*+ glueball masses in QCD3. All masses are in units of the square root of the 
string tension. 



State 


Lattice, N ^ 00 


Our prediction 


Difference, % 


2++ 


6.88 ±0.16 


6.72 


2.4 


2-+ 


6.89 ±0.21 


6.72 


2.5 


2++* 


8.62 ±0.38 


7.99 


7.6 


2-+* 


9.22 ±0.32 


7.99 


14 


2++** 


10.6 ±0.7 15 


9.26 


13 


2++*** 


— 


10.52 


— 



or 



(Tr (5VaV)^ Tr (5V5V) 



\x-y\- 



y' 27r44(:E - yY 



J2 (Af„M„)^/2e- 



(i\/„+M„)|x-y| 



(89) 



to obtain a paradoxical result: masses of 2++ states are given by the sums of two constituent 
masses M„, i.e are the same as masses of 0++ states! The resolution of this paradox lies 
apparently in the fact that we use local operators to probe extended objects. As a result. 



the correlator in (88 1 (or (81 1 ) describes not only spin-2 (or spin-0) resonances but states of 
other spins as well. Going back to our discussion of 0+^ sector, one may notice that in (82 1 



we identified as O^"*" only those states with one of the constituent masses equal to Mi (i.e. 
these states are of the form Mi ± M„, n = 1, 2, . . .). It should now be clear that 2++ states 
are described by the series with one of the constituent masses equal to M2, 4++ states - by 
series with one A/3, etc (we will say more on higher spin states in the next subsection). 
Thus our prediction for the masses of 2++ (as well as 2 ^) states is: 



M2++ 
M2++. 
M2++. 
M2++. 



= M2 + M2 = 8.42m 
= M2 ± M3 = 10.02m 
= M2±M4 = 11.61m 
= M2 + M5 = 13.19m. 



(90) 



In Table |5] this data is compared with the available results of lattice calculations. We see 
that like in 0+^ case we have excellent agreement for the lowest resonance while for excited 
states, the discrepancy is again at the 10 — 15% level. Let us note however that the lattice 
values for 2++* and 2 '"* states are quite different numerically (even though formally they 
are within the error bars) and this allows us to question their validity. This can be compared, 
for example, to the lattice predictions for 2++ and 2 ^ states which are almost identical 
(as they should be, due to parity doubling) and so seem to be reliable. Furthermore, as 
in the O'^''^ case, we notice that the lattice value for 2 ^* state is numerically much closer 
to our prediction for the 2++** state, while the lattice value for 2++** state is close to 
our prediction for the mass of the 2++*** state. Therefore, it is interesting to consider the 



^^Mass of 2++** state was computed on the lattice for 5(7(8) only I29| . The number quoted here was 
obtained by a simple rescaling of SU{8) result. 
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Table 6: Same as Table [5] but with lattice data relabelled as explained in text. 



State 


Lattice, N ^ oo 


Our prediction 


Difference, % 


2++ 


6.88 ±0.16 


6.72 


2.4 


2++* 


8.62 ±0.38 


7.99 


7.6 


2++** 


9.22 ±0.32 


9.26 


0.4 


2++*** 


10.6 ±0.7 


10.52 


0.8 



Table 7: 2^ glueball masses in QCD^. All masses are in units of the square root of the 
string tension. 



State 


Lattice, iV — > oo 


Our 


prediction 


Difference, % 


2+- 


8.04 ±0.50 




8.76 


8.6 


2— 


7.89 ±0.35 




8.76 


10.4 


2+-* 


9.97 ±0.91 




10.04 


0.7 


2—* 


9.46 ±0.66 




10.04 


5.6 



possibility that lattice states have not been properly identified and that proper labelling 
should be as follows: the lattice mass of the 2++* state can probably be considered as a 
reasonable approximation to the mass of a true 2++* state, but the lattice 2 ''* state should 
be identified as the true 2++** state. Finally, the lattice 2++** state should be compared 
to our 2++*** state. The comparison of our analytic predictions with corresponding lattice 
results after such relabelling is presented in Table |6] Once again, the numerical agreement 
is excellent. 

Let us proceed now to 2 states. For these states we may take Tr {dJd^Jd^ J) as the 
relevant probe operator. By computing the correlator of two such probe operators as we did 
before, we obtain that the masses of 2 resonances as sums of three M„'s, i.e. formally 
this is the same result as we have previously found for states. However by the same 
kind of reasoning as we had in the 2++ vs. 0++ case, we believe that the mass series for 
2 resonances should begin with Mi + M2 + M2. Thus our prediction for 2 sector of 
the theory is 



(91) 



For these states, there are not that many available lattice data. We list these together 
with our predictions in Table [7] Again we can see that agreement is reasonably good 
especially given the much larger error bars for the spin-2 lattice data as compared to their 
spin zero counterparts. 



2 — 


= Ml + M2 + M2 = 


= 10.99rn 


2--, 


= Ml + M2 ± M3 = 


= 12.59TO 


r 

2 — ** 


= Ml + M2 + M4 = 


= 14.18m 
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4.3 Higher spin states and Regge trajectories 

Unfortunately, at the moment, there are no lattice data available for glueballs with spins 
higher than 2 and it is our hope that future lattice simulations will address this issue. 
However, based on the spin-0 and spin-2 cases discussed above, it is possible to envisage 
natural higher spin generalizations of our results. One might expect that the quasi-Gaussian 
approximation to the vacuum wavefunctional begins to become insufficient to account for 
the masses of higher spin states. In the context of this paper, we have no way to test this, 
and simply provide the reader with the predictions. In particular, the masses of states with 
even spin and C — + are given by the sums of two constituent masses Mk according to the 
following rule {n = 0, 1, 2, . . .): 

• 0++ and corresponding resonances: Mg++. 

• 2++ and corresponding resonances: Af2++» 

• 4++ and corresponding resonances: M4++. 

• 6^+ and corresponding resonances: Mg++. 

or, in general 

Mj++,^ = Mj/2+1 + Mj/2+i+n, J = 0, 2, 4, 6, ... . (92) 

Here, for convenience, we choose to write these expressions for positive parity states only. 
One should remember that, by parity doubling, the masses of corresponding negative parity 
states are identical (unless J = 0) . 

Similarly, for the even spin resonances with C = — we may write: 

• and corresponding resonances: Mg--." = Mi 

• 2 and corresponding resonances: M2--." = Afi 

• 4 and corresponding resonances: M4--." = Mi 

• 6 and corresponding resonances: Afg--." = Mi 
and, in general 

Mj-,^ - Ml + Mj/2+i + Mj/2+i+n, J = 0,2,4,6, ... . (93) 

Once again we write these expressions for P = — states; the P — + results are identical. A 
few explicit examples of how Eqs. (92 93 1 work are given in Table [SJ 



= A/i 


+ Afi+„; 


= A.f2 


+ A/2+„; 


= Afg 


+ Af3+„; 


= Af4 


+ A'/4+„; 


J = 0,2,4,6 



ill + Afi+„; 


il2 + M2+„; 


lis + Af3+„; 


lU + M4+n; 



Given these generalized formulae for the masses of higher spin states, we can represent 
them graphically in the form of a Chew-Frautschi plot (J vs. M'^/a) and attempt to identify 
Regge trajectories. This is done in Figures 1 and 2. It should be noted, however, that to 
the best of our knowledge, there is no fundamental reason for the existence of linear Regge 
trajectories in the glueball sector of Yang-Mills theory. Nevertheless, our plots indicate that 
it is still possible to draw nearly linear Regge trajectories. 

We would like to conclude this section by pointing out that our discussion of the mass 
spectrum of 2+1 dimensional pure Yang-Mills theory is certainly not complete. In particular, 
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Table 8: Higher spin (4++, 4 ,6++, 6 ,8++) glueball masses in QCD^. All masses are 
in units of the square root of the string tension. 



State 


Our prediction 


4++ 

4++* 
4++** 


9.27 
10.54 
11.80 


6++ 

6++* 

6++** 


11.81 
13.07 
14.32 


8++ 

8++* 

8++** 


14.33 
15.59 
16.85 


4— 

4—* 


11.32 
12.59 


6— 
g— * 


13.86 
15.12 



current computational scheme based on the use of local probe operators is well suited for 

glueballs with PC = ++ or PC = . However for resonances with PC = — \- or PC = H — 

any probe operator with such quantum numbers is essentially nonlocal and it is not clear 
at the moment how our results generalize in that case. We were able to appeal to parity 
doubling phenomenon to circumvent this difficulty for states of nonzero spin J (although it 
is still interesting to compute masses of J^ and J ^ directly), but for ^ and 0"' we 
cannot do even that. It should be noted nevertheless, that it is still possible to pick certain 
combinations of constituent masses M„ to approximate existing lattice data on '" and 
0^ , but theoretical reasons behind this are not clear. 

Also we do not discuss here states of odd spin. However we would like to point out that 
existing lattice data for the lowest 1^+ states 



Mi++ = 9.98 ± 0.5 Vct 
Mi-+ = 10.06 ±0.4 y^^ 

can be approximated (up to 5%) by 

Mi++ = Ml + Ml + Ml + Ma = 9.50 y/^ 

while lattice values for 1*~ 

Mi+- == 9.43 ± 0.75 y^ 
Ml-- = 9.36 ± 0.60 V^ 

should presumably correspond to 

Mi+- = Ml + Ml + Ml + Ml + M2 = 11.55 y/a. 



(94) 



(95) 



(96) 



(97) 
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Figure 1: Chew-Frautschi plot of the large N glueball spectrum. Black boxes correspond 
to J"*""*" glueball resonances with even spins up to J = 12. Solid red lines represent linear 
Regge trajectories. 

As a possible explanation of such a large difference (about 20%) in the last case we may 
notice that in general it seems that C — — states are heavier than C — + states of the 
same spin (c/. 2++ and 2"' , for example). This is not the case with 1*~ lattice data and 
therefore we may question its validity. 

4.4 Fine Structure and Asymptotic Properties of Mass Spectrum 

As was extensively discussed in previous sections, glueball masses in 2 + 1-dimensional pure 
Yang-Mills theory are given by expressions involving zeros 72, i of Bessel function J2{x). For 
example, equation (92) describing J"^"*"* (for even J) resonances can be rewritten in terms 
of 72,1 as 

Mj+ + ,r, = — (72, J/2+1 + 72,,7/2+1+n) • (98) 

We can use now the following asymptotic expressions for the zeros of Bessel functions 



m 



lu,i 



1 1 

7r| n+-v- - 



4z/2 



?7r(n+ \v- \) 



O 



l2- 



TT I n + 



15 



87r(n+f) 



O 



(99) 



(100) 



to obtain asymptotic expressions for highly excited J++* states (with even J and large n) 
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Figure 2: Same data as in Figure 1. Red lines represent polynomial data fit with second 
order polynomial J = ao + ai I ^ ) + a2 ( ^ ) • Note also the different way to connect 
data points compared to Figure 1. 



Mo4 



— m— n + 



7 
4 



72,1 



o 



M. 



2++* 



11 



n 
m— \ n , 
2 \ 4 



72,2 
TT 



o 



,, TT / 2J+7 72„//2+l 



O 



(101) 



We can see now that for large excitation level n, the glueball mass spectrum with given J^'-^ 
quantum numbers becomes equidistant with mass splittings 



AM = Afj_^^.(„+i) - Mj++.^ = m-. 



(102) 



For large values of spin J we can further simplify Eq.(101l if we replace 72,7/2+1 by its 
asymptotic value. We get 



M,..,.=mUn + J+l)+o[K^j 



(103) 
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Figure 3: "Band" structure of the mass spectrum. 



This is an interesting expression since it tells us that there is approximate mass degeneracy 
in the theory. By this we mean that the mass of a state with given large values of n and 
J is approximately the same as masses of states with the corresponding excitation number 
and spin equal to n + i and J — z respectively {i — ±2, ±4, . . .). Actually, since Bessel zeros 
approach their asymptotic values very fast we can see such approximate degeneracy even 



for states of low spin J. For example, as was found in (82 1 and (90) masses of 2++ and 

0++** states 

Mo++-. = Ml + Mg = 8.38m 

M2++ = M2 + M2 = 8.42m 
are very close numerically, the difference being about .047%. Similarly for 



(104) 



Mo++- 
M2++. 



A/i 
M2 



A/4 = 9.97m 
M3 = 10.02m 



(105) 



the difference is about 0.52%. 

We thus see that the mass spectrum of the theory has a "fine structure" with numbers 
of minutely separated mass eigenstates grouped together into "bands". It is interesting 
nevertheless to investigate their structure a bit further. In general, 2n's "band" is made out 
of the states with the following mass content: Mn + Af„, Af„„i + Mn+i, ■ ■ ., Mi + A'/2„-i 
(and similarly for "bands" with odd number 2n + 1). To see their general features, it is 
convenient to work out an explicit example; say, for the 50th "band" we have 



Ml + A/49 = 80.70893m 
Af2 + A/48 = 80.77882m 

A//23 + Af27 = 80.872686m 
M24 + Ar26 = 80.872792m 
M25 + Ar25 = 80.872827m. 



(106) 
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From this example it is not difficult to see that the heaviest state within a "band" is the 
state with the most equal masses (M25 + M25 in this example), while the lightest state 
is the one with the most different constituent masses (Mi + M49). The "bands" do not 
overlap: say, the most massive state in the 49th "band" is M24 + M25 — 79.301557TT- and so 
the distance between two "bands" (about 1.5m) is about 10 times the width of the "band" 
(about 0.16m). Finally, as can be seen from ( |106 l, mass levels are very dense near the top 
of the "band" and spread further apart towards the bottom. Using Eq.(100l it is possible 
to derive the needed asymptotic expression for the mass splittings within a "band" . The 
highest level in the 2n'th "band" is 



M„ + M„ 



TO 
'2 



27r n 



30 



8^(n+|) 



O 



1 



(107) 



Next levels are 



M, 



n-^a 



M„ 



m, 



2n{n+ - 



30 



30 



O 



8n{n+l) 8^(n+f) 
and therefore the distances between these levels and the highest level in the "band" are 

15to, 



(108) 



AM„ 



'Sn{n+lf 



a 



O 



(109) 



for a = 1,2,3.... Figure 3 gives graphical representation of the "band" structure. 

A similar analysis can be done for states with other J^'-^ quantum numbers. For brevity, 
we present here only the asymptotic expressions for J state masses: 



Afn 



to — 
2 



, 7 72,1 +72,1 

4 TT 



n 



Mn 



11 72a + 72.2 

4 TT 



O 



,. '^ /^ , 2J + 7 72,i+72,j/2+i 



and for large values of spin J equation ( 110 1 further simplifies to 

1 1 



M 



J- 



^ I , 7 , '^ , '^2,1 



o 



n' J 



O 



(110) 



(HI) 



It is interesting to note that asymptotic expressions (101 1 and (110 1 are reminiscent of 
empirical observations made in |28j . 
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Also, if wc include all of the spin states, the general structure appears to be reminiscent 
of the spectrum of a string theory. Although degeneracies are not exact, the bands that we 
have discussed here appear to be identifiable with the levels of a string spectrum, and there 
is essentially an exponentially rising density of states. This is the sense in which this theory 
gives a Hagedorn density of states. The fact that degeneracies are not exact, particularly 
at the lowest mass levels, indicates that this is not a free string theory. 

5 Conclusion and Open Questions 

Based on the formalism of "corner variables" [H |6] we have presented an evaluation of 
a quasi-Gaussian wave functional for Yang-Mills theory in 2 + 1 dimensions containing a 
very non-trivial kernel compatible with an infinite number of relativistic point-like resonant 
states. The masses of these states are all given as combinations of Bessel function zeros. The 
actual numerical values thus obtained are in remarkable agreement with lattice simulations 
[8]. The resulting picture of glueballs is based on gauge invariant constituent states, and is 
rather of an open-string than closed-string (Wilson loop) like nature. There are a number 
of outstanding technical issues, and we have been careful to point out throughout the paper 
the status of each step. 

The perceived success is based heavily on comparisons to the existing lattice data, the 
most reliable of which are for the lowest lying low spin states. As one goes to higher spin 
states, there is certainly a lot of physical phenomena that could cloud the picture, such 
as the mixing of states. Although we have presented some results for non-zero spin, we 
should perhaps remain cautious, as there are little data to compare to. It is also possible 
that for some observables, the quasi-Gaussian form of the vacuum wavefunctional becomes 
insufficient in some way. It would be useful, for example, to construct the actual glueball 
wavefunctionals and compare their energies to the masses that we have derived here. 

In conclusion, we believe that our results demonstrate the importance of "corner vari- 
ables" jH |B] in the pure Yang-Mills sector. As mentioned in the introduction to this paper, 
these variables are not only useful in 2 -I- 1 but also in 3 -I- 1 dimensions [7j. Many other 
purely theoretical questions remain: understanding of non-perturbative renormalization, 
better understanding of the constituent picture, exploration of the lattice formulation in 
terms of "corner variables" , understanding of a manifestly covariant formulation etc. Obvi- 
ously on a more pragmatic level, it is important to understand the inclusion of quarks (and 
the emergence of mesons and baryons) in this approach in both 2-1-1 and 3-1-1 dimensions. 
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6 Appendix A: Regulated Computations 

In this Appendix we want to discuss the eigenstates of the kinetic operator 

In particular we are interested in the action of T on a special class of local operators 
On — /9J(A)"9J. Co = J dJ°'dJ'^ has been considered previously in detail in U and 
we wish to generalize to arbitrary n. Since we will be unable to do this infinite number 
of calculations, we will focus on the general structure of TOn- Ultimately, we will find a 
number of issues in these calculations that we are unable to resolve. The issues include 
the use of the holomorphic invariant point-splitting regulator and precisely how to remove 
the regulator, as well as renormalization and normal-ordering issues F^ The reader should 
interpret this Appendix as a survey of the issues involved, as well as motivations for eq. 



(56 1; this equation is not derived. 

We start with a very important technical point. The action of the operator 



on any local operator-valued function of J is always calculationally straightforward. Ti 
simply acts as a number operator counting J's. For example, the action of Ti on Oq = 
J dJ"'dJ'^ is trivial and we simply get 

TiOo = 277100. (114) 

More generally though, the action of Ti on ©„ will not result in a holomorphic invariant. 
On the other hand the action of the other part of the kinetic operator involving 

T,.™£n.'(.,.)^j^ (115) 

is very subtle. For Oq we have 



^^Another route to this computation is the one followed by Karabali, Kim and Nair in their computation 
of the action of the kinetic operator on Co- Their computation is based on an explicit insertion of Wilson 
lines in the point-split form of our Oq operator. This may be taken to mean a different basis of operators. 
There is apparently a mismatch between our computation and theirs, which can be traced to the ultralocal 
form of our On operators. The Karabali-Kim-Nair type of computation has not been generalized in this 
paper. Obviously it is important to pursue that approach to check the results presented in this paper. 
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and so 



T2O0 = -2m / 6"'' [dln''''{z,w)]^ 



(117) 



To compute the right-hand side of Eq.(117), we use the point-spht expression for fl {z,w) 

a 



n-''{z,w)^D'::Ara{w,z) 



(118) 



where D^ is the adjoint representation of the holomorphic-covariant derivative, and [4] 



Ara{w,z) = 



tt{z — w) 



[Sra - {Hiw)H-\z,w))^^e-^ + 



(119) 



with a = "' , and e ^ 0. Now one easily obtains (summation over repeated index is 



understood) 
and finally 



[d',n^^{z,w)]^ 



T2O0 



Cyi {dim G) 
CA {dim,G) 



(120) 



(121) 



TT 27re2 

From this we see that the only effect of T2 acting on Oq is to give a divergent normal-ordering 
correction. We may now define 



:Oo: 



drdr 



CA (dimG) 
TT 47re^ 



and for normal-ordered operator : Oq : we get 

T : Co : = 2to : Co : 



(122) 



(123) 



In what follows we will not keep track of such divergent terms explicitly, however it should 
be understood that the rest of 0„ operators have to be normal-ordered in a similar way. 
Let us proceed now to the operator 



Oi = / dJ''iA)''''dj'' 
The action of Ti on this is elementary 



Qjaj^abg2jb_ 



Ti Oi = 2m / dJ^id) dV" + 3m / dJ" ( -i — f'""'J'' ] d'j" 



(124) 



(125) 



Note however that this result is not invariant under holomorphic transformations. What 
we want to show now is that when T2 acts on the 3 J-part of Oi (the action on the 2 J-part 
is trivial and gives only a normal-ordering constant as was discussed above) it generates an 
extra term with two J's which restores holomorphic invariance, i.e. 



and therefore 



TOi ^{Ti+T2)Oi = 3mOi. 



(126) 
(127) 
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But before we proceed to show that ( 126 1 is correct we want to make a technical remark 



which will significantly simplify calculations. As can be seen from expressions (1181 and 
( 119 1 every time dz acts on il°-^{z, w) it actually acts on the exponent in A^a and thus pulls 



out a factor 1/e. Therefore expressions hke [^^^^''''(z, w)] _ , [dzdyjn°-''{z,w)] _ , eic will 
lead to divergent normal-ordering terms only (c/. Eq.( 120| , for example). Such expres- 
sions will appear when S/SJ°'{z) partial derivative acts on J's with leading antiholomorphic 
derivatives, like dJ or d^J in the case of Oi operator. Therefore, as long as we are not 
interested in explicit form of the normal-ordering terms, we can consider action of S/SJ°'{z) 
on "bare" J's only. Therefore 



dJ°-[z) \ CA 



and 



(501 
6J<'{z)5J''{w) ^ y \a 
from which we immediately obtain 

T2 d === m 



^fae\ f^Q2je gj(^^^ ^ w) - 9Jf d^J{z - w)) 



(128) 
(129) 



^-f-'^j [B^r [d^^n'^\z, w)] ^^^ - Br [Bln'^\z, w)] ^^J . (130) 

To proceed further we need to evaluate [Byjn°'^{z,w)'\ _ and [9^ri''''(2;, tu)! _ . A 
straightforward computation gives (we keep only finite terms which do not depend on reg- 
ularization parameter e) 



b/., „..\^ _ ^ a/T-,T\ba_ 1 dBj'"^ — ^ "" ^' ^ ^^''' 



[B.^n^\z,n.)]^^,^^-B{Djr^^^ 



and 



[dln^'iz^w)]^^^ 



Ta d = m 



ZTT 27r 



2 jba 



2tt ca 

1 TT 

27r CA 



B{jjr 
bvj)'" 



_■ jbae \ fQ2je ggjba _ Qje gg2jba\ 

2CA ' ^ 



(131) 

(132) 
(133) 



Finally, we may use the identity f<>°-'^j'"^ = icA J'^ to verify that (1331 is indeed equivalent 



to (126) 



To summarize, what we have shown is that the action of kinetic energy operator on Oi 
gives a gauge and holomorphic invariant result 



r : Oi := 3m : Oi : . 
Finally, let us consider one further example 

It is convenient to rewrite this as 

©2= ({Dd^jYiDB^jf. 



(134) 



(135) 



(136) 
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Proceeding now the same way as we did for Oi, we obtain 



and 



SJ^iz) CA 



, /,^\, , =2—1 i-if-' ~dH{z - w)) (D B^jY + 

CA 



(137) 



(138) 



Now we may write for the action of T2 on O2 



T2O2 =2m- \ / {-if'^' [dl^^\z,w)]^_^J) {Dd'j) + 



(139) 



J CA 



Since [9^^„f2''''(z, w)] _ was computed previously in (1321, we immediately obtain for 
the first term in ( 139 1 



2m- / {~^r^'' [dln^\z,w)]^_^J {Dd^J)l = mj {dd'jy {D B^J)' 



Similarly for the last term in (1391 we need expression for [il°^(2;, i«)] _ which is 

1 



[n'^\z,w)] 



2ti^ ' 27r 



ba 1 Qj,. 1 TT 



27r CA 



{33) 



ha 



(140) 



(141) 



and therefore 
2 



2m 



TT 



CA 



-if"-" B^3^) [n''''{z, w) 



if" diJl) = -m 



1 / TT 



TT \CA 



{dl3B^3D3Y^ 

(142) 
where summation over the repeated index a on the right-hand side of this expression is 
implied (this essentially amounts to taking trace in adjoint representation). Finally, for the 



second term in (139) we will need 



[D'::^dln'^\z,w)] 

from which we obtain 



= L\dB''D3--B^D-'3 
27r I 3 



da 



(143) 



2m^y {~^f^-S'3^ {[D^!^ Bln<^\z , w)] .=J = -^ 



d^3^Dd^D3-^a^D^3 



(144) 
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Now if we collect all the results in equations (140), (142) and (144) together and after some 



simple but somewhat lengthy transformations we obtain the final expression for the action 
of T2 operator on O2 state 

T,02^jUJ {dd'jy {dd^jy - (145) 

-[j,9V]yf-[j,a2j] 

CA J \CA 

— {[dJ,d^J])\Ddjf\ 
CA J 

which combined with the straightforward result for T1O2 

Ti 02 = "I |2 /" {dd'^jy {dd^jy - (146) 



CA J \CA 



finally gives 



T02^{A^^^m02+'^^^j m~d'j]y{Ddjy. (147) 

As a consistency check one may notice that we obtain a gauge and holomorphically invariant 
expression as expected. 

A few comments about this result are in order, as it is not of the form that we have 
expected. First of all we see that an extra gauge invariant operator J ([9J, 9^J]) (D9J) 
mixes with 02- Such operator mixing is certainly expected to be the general feature of the 
theory and therefore we may write 

TOn - E,,On + . . . (148) 

where ". . ." indicate operators which mix with On under the action of T. It should be noted, 
however, that such operators, even though they have the same mass dimension as ©„, are 



of higher order in dJ. By this we mean that only 0„ in equation ( 148 1 contains terms 
quadratic in dJ while the rest of the operators have at least three 9J's. This means, in 
particular, that such operator mixing has no influence on derivation of Schrodinger equation 



presented in this paper and therefore detailed knowledge of ". . ." terms in equation (148) is 
not really necessary for our purposes. 



More worrisome is the numerical coefficient of O2 in ( 147 ) . As can be seen from ( 147 1 , 
we may think of O2 as an eigenstate (up to an extra operator that mixes in) of kinetic 
energy with eigenvalue (4 — 1/3)to. 

We believe, however, that this result is not complete and the correct spectrum for the 
kinetic energy should be equidistant 

TOn = {n + 2)mOn + ■■■■ (149) 
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The reason why we have not obtained the expected answer in the case of O2 (and presumably 
for On>2) niay be that we have not taken proper account of the nonlocal character of the 
theory. Thus the calculation that we have outlined above needs to be reconsidered carefully. 
We do not at this time have a consistent understanding of all of these issues as the required 
calculations are tedious, but we present remarks below. In the text of the paper, we use eq. 



( 149 1 without further apologies. 



6.1 Further discussion: regulator/renormalization issues 

In deriving results for Oq, Oi and O2 operators we have been using expression for Q.°'^{z, w) 



with Kra{w,z) given in equation (1191. This expression for Ara(w,z) is in fact a simplified 
form of the exact regulator j4] 

Ka{w,z) = -^ ^ [5ra-{H{w)H~\u,w)H{u)H-\z,u))e-^+... (150) 



where u = ^Az + w) and u = i(z + w). It is easy to see that (1191 can be obtained from 



this expression if we put z — w 'vn (^H{w)H ^{u,'w)H{u)H ^{z~uyj . Therefore (1191 is 



equivalent to ( 150 1 up to the 0{z — w) terms. Part of the problem uncovered in the previous 
section then may have to do with subtleties involving the assumed form of A^a- It seems 
that for the study of the action of T2 on local operators, such simplification is not justified 



and exact expression ( 150 1 should be used. Although it becomes extremely difficult to do 
calculations with this Ara{w, z), we have considered this effect on Oq and Oi. There appear 
to be additional contributions (of the form e/e) even to their eigenvalues 

TOo-^2(l-^^mOo (151) 

What does this effect illustrate? There are certainly very subtle effects coming from the 
point-splitting regulator. We have not been able to arrive at a consistent method for dealing 
with these effects. In addition though, we should note that these are all formal expressions 
- we have not carefully defined the operators or the coupling, particularly in the infrared. 
It is not clear to us that such a procedure is known in this or any other (non-trivial) theory. 
Let us be a little more explicit about the issue raised here. In performing computations 
for Oq, Oi and O2 we have kept track of e-independent terms only. In general, however, we 
may write 

T Oo - . . . + Too Oo + Tio e Oi + T20 e^ O2 + ■ ■ ■ 
TOi = ... + Tn Oi + T2ie02 + T31 e^ O3 + . . . 
T02 = ... + T2202 + T^2(0^ + Ti2e^O^ + ... ^ ' 



or in matrix notation 

TO,^T,,e^-'0,. (153) 

This is a schematic expression based solely on dimensional analysis (remember that [On] — 
(mass) "). It does not include, for example, the effects of mixing with operators other 
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than On- Tij here is an infinite dimensional number matrix which can be thought of as a 
representation of kinetic energy operator in the basis provided by operators 0„ . Evaluating 
the spectrum of T would then certainly amount to diagonalizing Tij. In some sense what 
we have done above is have found diagonal elements T„„ of this matrix and associated them 



with the spectrum of kinetic operator. In other words if one takes naive e — > limit in ( 152 1 
it looks like all terms with positive powers of e go to zero, Tij matrix becomes triangular 
and, consequently, diagonal elements give the spectrum. The proper procedure, however, 
would be to diagonalize Tij first and take e — > limit afterwards. Formally we may write 
the n-th eigenstate as a linear combination of an infinite number of local operators ©„ 

Q{n) ^ ^(n) g»-n q_ (^54) 

with some unknown numeric coefficients c^ . The eigenvalue equation 

rg(«) ^^(«)g(") (155) 

becomes now the infinite-dimensional matrix equation 

Tj, c|") = ^(") 4") (156) 

from which in principle we could find the spectrum. Of course, to do this practically would 
be a very difficult task. Note, however, that the e-dependence completely drops out from 
this expression. This means, in particular, that the spectrum of T does not depend on 



the regularization parameter e and we may take the formal limit e — *■ in (1541 after 
diagonalization. At any finite value of e the true eigenstates Q^"-* of T are certainly nonlocal, 
however we expect that they will collapse into local states as e — + 0, i.e. 

Qin)'^.,0„:. (157) 

Therefore for any practical purposes, like writing Schrodinger equation for the full Hamilto- 
nian Ti, we may think of the eigenstates of kinetic energy in terms of local operators : On '■ 
and we may now really write 

T :©„ :'="£;(") : ©„:+... (158) 

It should be noted however that this expression, even though it looks completely local, knows 
in fact about nonlocal character of the theory. This nonlocal information is contained in 
the spectrum E^"' which should be obtained by solving eigenvalue equations with nonlocal 



ansatz (154 1. 



To conclude this discussion, we would like to consider a simple example which further 
illustrates the intuitive arguments presented above. In [1] the following nonlocal generaliza- 
tion of Co was proposed 

2^°^=/ —dr{x)[H{x,y)H-Hy)T''dr{x). (159) 

Jx.y 

This state is essentially made of two 9J's connected by a straight Wilson line. It is not 
difficult now to expand this in terms of local operators with the result given by 

Q(°) = 7^o + eUi + ^e27^2 + j/U^ + ... (160) 

36 



where 



Tlr. 



ah 



dJ" D''d") dJ' 



(161) 



These ??.„ operators are not quite the same as On (except 7?.o = Co and TZi = Oi), however 
they coincide with the latter to quadratic order in J. Therefore from this e xpansion we see 



)(o) ^o 



that Qe is of the same type as we need to diagonahze T according to ( 154 1. Of course, this 



simple operator is not an eigenstate of T, however we may see the other essential ingredient 
of the above analysis, namely the fact that in the continuum e ^ limit, Qe collapses into 
local operator On- 



7 Appendix B: Vacuum Energy 

Here we discuss the vacuum energy. For this purpose, we need to evaluate the Schrodinger 
equation in the form 

T^^^AT^o = (^0 + • • O^o (162) 

where Eq consists of terms that survive as J" — » 0. Since the Hamiltonian involves at most 
two derivatives with respect to J, we can drop any J-dependence of the kernel for this 
computation and work with 



^0 — exp 



We then compute 



and 



SJ^iz) 



*n 



2cA'mr 



CAm' 



drK[dd/m'^]dr 



{dK[dd/m'^]dJ)l^o 



*n 



and so 



5Jb{w)5J'^{z) ' CATu 



m / r2"''(z,u)) 



-^ {K[dd/m^]d^)j'^^\z-w)S'''''<5Jo + .. 



6^ 



*o 



SJ''{w)6J''{z) 
We then deduce that the vacuum energy density is 



£o 



*0- 



Si 







6,,{K[dd/m']d')^n^'{z,w)l^ 



(163) 

(164) 
(165) 
(166) 

(167) 



CAITI 

To proceed, we write out the regulated fi'^'' keeping only those terms that are present in the 
J° ^ limit. This is particularly simple, and we find 

, r m) = .^^( n. Mil) rU"" = -^^ A"" r). 

eg\ 



n^l^{z,w)=''-^{D^,K{w,z)f-=''-^5-''d.^ . , 

TT TT 7r(Z — W) 



1-e" 



-|z-ujr/2e 



(168) 



In the expression for the vacuum energy, there is a leading 9^; performing one of these 
derivatives on fi^g^, we find 



fo 



771 ^ '■ i 'z 27re 



= —m{dimG) K 



dd 



dd 



1 



^-\z~w\'/2. 



m? 1 . 27re 



(169) 
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This expression is horribly divergent. Let us proceed however by considering the expansion 
for K 

OO 

K[x\x^^Cr,x''+\ (170) 

We then find 



En = - 



{dim G) m 



27re .^ "^ ' \ 2em^ 



n+1 






Note that the appearance of divergences to arbitrary order here conies directly from the 
fact that the vacuum wave-functional contains operators of arbitrarily high dimension. This 
basic fact plagues all such computations. In the case of the vacuum energy, it is possible 
to resum this series and extract sensible results (in particular, the asymptotic behavior is 
consistent with the UV perturbation theory). The basic observation is that 



(n+l)!t-2-n= / 
Jo 



dx a;"+^e-*^ (172) 



and so we have 



,3 °° 



£0 = (rf»^G)m y^^(_,). Hdxx-^h 



n=0 
-.3 /-OO 

dx xK(—x)e~ 



(dimG) m , .^/ x -tx 



(173) 
/o 

Thus the vacuum energy may be written as the Laplace transform of xK{~x). Removing 
the regulator corresponds to taken the asymptotic t = 2em? — > limit. 

7.1 Evaluation of Sq 

To proceed further, we need to supply information about K. We will consider several 



examples, culminating with the kernel (62 1 



7.1.1 K=l 

The simplest choice is to consider K = 1, which corresponds to replacing the true K by its 
IR limit. In this case, the Laplace transform is elementary, and we obtain 

which is quartically divergent (in powers of a momentum cutoff). 

7.1.2 Massless Boson 

If we consider the wavefunctional corresponding to a free massless boson, we have K{—x^ = 
\/ ^/x. Again, the Laplace transform is elementary, and we obtain 

pparticle dlTTl Lz , , 



This is also the UV limit of the true K; we see here that the scale m drops out of this result, 
as it should. 
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7.1.3 Bessel 



Finally, let us consider the kernel (62 1. K{—x) can be written as 

1 /2(4\/^) 



K{~x) 



/^/i(4y^) 



and so we need to compute 






(176) 



(177) 



Because we are interested in the small t limit, it is sufficient to consider the asymptotic 
a; ^ oo behavior of the integrand 



'/i(4VS) 



3 1 

128^^1 



0{l/x). 



We thus obtain 



YM _ dim G 
° ~ 4v^e3/2 



1 



3?Ti-v/e 3m e 



2v/2^ 
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+ 0{m^log{em)). 



(178) 



(179) 



Note that the leading divergence is the same as the UV result, while the subleading terms 
depend on the detailed form of the kernel. 

8 Appendix C: Vacuum Solutions 



Assuming the validity of the Riccati equation ( 59 ) 



K-~[K{L)]+LK' + 1 = 



(180) 



we now consider its general solution. This is a non-linear differential equation but it is easily 
solved by substituting K = —y'/2y to obtain 



Ly" + 2y' + iy = 0. 
Then, writing y = f{x)/^/L with x = 4\/Z leads to the Bessel equation 

nx)+nx)/x+{i-i/x')f{x) = o 

with general solution 

2/ - ci 4^ Ji (4 VZ) + C2 4^ Fi (4\/Z) . 
\/ L \/L 

Using the standard recursive formulae (with similar formulae for Y,^) 

J^-i{u) + J^+i{u) = ^Mu) 

J^-i{u) - J^+i{u) = 2J^(u) 



(181) 

(182) 
(183) 

(184) 
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the general solution of ( 180 1 takes the form 



i.(L) = ^^^Ml^^^±^l(l^. (185) 

VlCJi{4Vl)+Yi{4Vl) 



Here we have a one-parameter family of solutions, parameterized by C. The solution quoted 



in the text (62 1 corresponds to C — > oo. To see why we should take this value, we study the 
asymptotics of the solutions, paying close attention to the normalizability of the solution. 
It is particularly satisfying that the only normalizable solution corresponds to the confining 
vacuum. 

To study the asympotics, we note that L = dd/m"^ — > —k'^ /Am? is negative (choice of 
square root doesn't matter), and so we need the identities 

Ji{iz) = ili{z), J2{iz) = -hiz), 

Yi{iz) = -h{z) + -Ki{z), Y2{iz) = -il2{z) + -K2{z) (186) 

TT TT 

giving 



, , ,, 1 C'l2{AJ\L\)-K2{AJ\L\) , , 

V^C'/i(4v^) + ifi(4v^) 



where C = 7r(C + «)/2. We will consider real C" in what follows (this corresponds to 
real wavefunctionals) , but this is not a crucial simplification: the final conclusion will be 
unchanged. 

Now, asymptotically \L\ — )■ oo, this function behaves differently, depending on whether 
C" = or not, because only for that value of C" will the K^s play a role: 



1=, /or|C"|7^0 



l\L\ 

For small \L\ on the other hand, the Kn& will dominate, unless C" -^ ±oo. Thus we find 

I 1, JOT C — ±CX3. 

Thus there are three types of solutions, in terms of their asymptotics. Let us discuss them 
in turn: 

1. C" = 0: for this solution K{—\L\) < for each value of \L\ meaning that this solution 
is completely non-normalizable (modes of all momenta are non-normalizable) ; 

2. C" ^ 0, C" 7^ ±00: for these solutions ii'(— |L|) changes sign from negative to positive 
at some finite value of |L| meaning that it is normalizable in the ultraviolet but non- 
normalizable in the infrared; 

3. C" — ±cxd: for this solution iir( — |i|) is everywhere positive meaning that this solution 
is completely normalizable. 
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To summarize, as one can see from this discussion the only sokition of Riccati equation 



( 180 1 which is normahzable in both UV and IR corresponds to C" — ±00 and is given by 



K{L) = 



1 J2(4\/L) 



/ZJi(4\/L)' 

9 Appendix D: Bound states 



(190) 



For a free scalar, the Hamiltonian is 

1 



Us 



im' - V^ 



and the vacuum state is 



*o = cxp - 



1 



4>{x) 



V fc^ + m? 



<i>{v) 



The kernel is therefore Ks{k) — y k'^ + m? and the inverse kernel is 

1 



K7'{k) 



yp 



In 2+1 dimensions, we can also deduce, for spatial separation 

D{x - y) = {4>{x)4>{y)) 



m\x—y\ 



(191) 



(192) 



(193) 



(194) 



27r|f -y| 
which is equivalent to a propagator with single pole. This functional form is what we found 



in the factorization of the O'''"'" correlator (81 1 



In our case, we can also see this physics in momentum space. Now, in our case the 
inverse kernel takes the form 



K-\k) = 



1 



(195) 



fc2 + Af 2 

Compute now the Fourier transform of {K^^{x — y)Y ■ We claimed in the paper that it 



had the form of (194) 



{K-'Y[k] 



(27r)2 



K-'^{p)K-^{k~p). 



Given the expansion of K ^(fc) we may rewrite this as 

dPp 1 



{K-y[k]= ...+ 



{2t:Y p^ + Mf {k-p)^ + Alf 



(196) 



(197) 



so essentially we have to compute a simple two-dimensional scalar loop integral. This is a 
simple computation and, for example, for equal constituent masses i = j we get 



(K-'m 



1 



ArcTan 



fc/2 



v/fc2 + Ulf 



2TTk 



(198) 
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Now we immediately see that this expression describes a physical single-particle state of 
mass M = 2Mi: the leading square root factor is the same as for a free scalar, as above, 
while the multiplicative function of momentum does not seem to have any unwanted analytic 
properties. 
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